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ABSTRACT 

The  propagation  of  the  electromagnetic  waves  in  vacuum  is 
discussed  in  terms  of  the  Lorentz  Transformation.  A  simple,  but 
rigorous,  derivation  of  the  Lorentz  Transformation  is  given. 

Then,  the  observed  doppler  shift,  seen  by  a  moving  observer,  of  a 
moving  target  is  derived  in  terms  of  inertial  position  and 
velocity.  This  radar  doppler  equation,  given  without 
approximation,  is  suitable  for  analysis  of  doppler  data  for 
satellite  orbit  determination. 


I.  INTRODUCTION 


Propagation  of  Electromagnetic  Wave  is  rigorously  described 
by  Maxwell's  equations.  The  description  of  Electromagnetic  Wave 
propagation  in  one  reference  frame  is  related  to  the  description 
of  the  same  wave  in  a  second  reference  frame  by  the  General 
Lorentz  Transformation  (GLT),  which  is  a  direct  consequence  of 
the  Special  Theory  of  Relativity  (TSR).  The  GLT  encompasses  both 
a  change  in  the  bases  of  the  reference  frames  and  a  relative 
velocity  of  the  two  reference  frames.  The  GLT  is  a  unitary 
transformation  of  a  four-vector,  i.e.,  preserves  the  length  of 
the  vector.  The  GLT  has  all  the  usual  properties,  suitably 
generalized,  of  unitary  transformations. 

In  the  following  development  we  will  discuss  the  basic  ideas 
that  lead  to  the  GLT,  in  a  superficial  but  rigorous  way.  That 
is,  the  results  are  correct,  but  many  of  the  fundamental 
consequences  will  not  be  addressed.  These  broader  consequences 
of  the  GLT  are  addressed  in  many  classical  texts,  and  need  not  be 
repeated  here.  However,  though  the  final  results  given  here 
cannot  be  considered  new,  they  are  not  explicitly  given  in  any 
text.  Therefore  we  hope  to  provide  a  complete  and  rigorous 
result,  with  sufficient  detail  that  it  can  be  applied  directly  in 
practical  application  and  without  ambiguity. 

The  radar  Doppler  equation  involves  four  reference  frames: 


1.  The  frame  of  the  transmitter  at  the  time  of 
transmission,  in  which  the  transmission  frequency  is 
observed . 

2.  The  frame  of  the  target  (satellite)  at  the  time  of 
reflection  or  retransmission  by  a  transponder,  the  same 
thing  in  the  context  of  this  discussion. 

3.  The  frame  of  the  receiver  at  the  time  of  reception,  in 
which  the  received  frequency  is  observed. 

4.  The  geocentric  frame  in  which  the  target's  trajectory 
(ephemeris)  is  described. 

The  Lorentz  transformation  is  a  mathematical  description  of 
the  t^o  basic  tenents  of  the  Theory  of  Special  Relativity  (TSR), 
as  stated  by  Einstein  in  1905: 

1.  All  the  laws  of  nature  are  identical  in  all  inertial 
systems  of  reference.  An  inertial  reference  system  is 
one  in  which  a  body,  not  acted  upon  by  external  forces, 
proceeds  with  constant  velocity. 

2.  The  velocity  of  propagation  of  interaction  between 
particles,  which  is  the  same  as  the  velocity  of  light 
in  empty  space,  is  the  same  in  all  inertial  reference 
frames . 

The  TSR  denies  the  existence  of  any  fundamental  reference  frame. 
Therefore,  none  of  these  reference  frames  is  any  more  fundamental 
than  another.  The  principal  objective  of  investigating  the  Radar 
Doppler  Equation  is  for  determination  of  the  target's  trajectory 
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(orbit  determination)  based  on  the  observed  Doppler  shift. 
Therefore,  for  the  following  discussion  we  will  consider  the 
geocentric  frame  as  the  basic  frame,  and  seek  to  express  the 
observed  doppler  shift  as  a  function  of  the  geocentric  position 
and  velocity. 

II.  THE  LORENTZ  TRANSFORMATION 

We  give  here  a  simplified  but  complete  derivation  of  the 
Lorentz  transformation.  There  are  a  number  of  lucid  developments 
of  the  Lorentz  transformation  [Refs.  1,  2,  and  3)  that  reveal  far 
more  general  properties.  However  for  our  purposes,  precision 
metrology,  what  follows  here  allows  us  to  be  precisely  clear 
about  each  quantity,  when  we  come  to  practical  calculation. 

A  reference  frame  is  understood  to  consist  of  a  system  of 
meter  sticks  and  clocks  to  define  coordinates  (x,y,z)  and  time 
(t).  Now  in  a  reference  frame  we  can  imagine  relating  the  space 
and  time  coordinates  by  using  the  travel  time  (t2  -  tj)  of 
electromagnetic  propagation  between  two  points  (x2  ,  y2  ,  z2 )  and 
(*l»  Yi .  zi)  as: 

c2(t2  -to2  -  (x2  -  *i)2  -  (y2  -  yi)2  -  (z2  -  zi)2  -  o  (i) 

A  fundamental  concept  in  the  TSR  is  an  event  defined  by  the  four 
quantities  (t,x,y,z).  Now  we  define  the  (proper)  distance,  i.e., 
the  interval,  between  two  events  (s)  as: 


(2) 


s2  =  c 2 ( 1 2  -  to2  -  (x 2  -  xt)2  -  (y2  -  yi)2  -  (Z2  -  Zi)2 

or  as  differentials: 

ds2  =  c2dt2  -  dx2  -  dy2  -  dz2  .  (3) 

The  electromagnetic  propagation  in  a  vacuum  is  said  to  have  an 
interval  s  =  0  as  given  by  (1).  It  follows  from  the  second 
tenant  of  the  TSR  that  if  the  interval  between  two  events  is  zero 
in  one  inertial  reference  frame,  then  it  is  zero  in  all  other 
inertial  reference  frames. 

To  make  this  more  concrete,  consider  an  inertial  reference 
frame  (K')  moving  with  uniform  velocity  (V)  with  respect  to  an 
inertial  reference  frame  (K)  along  the  x  axis  as  illustrated  in 
Fig,  1.  The  differential  interval  between  two  events  (ds)  in  the 
K  frame  is  given  by  (3).  The  same  two  events  viewed  in  the  K' 
frame  would  have  an  interval  (ds  ' )  : 

(ds')2  =  c2 (dt ' ) 2  -  (dx ' ) 2  -  (dy')2  -  (dz')2  .  (4) 

We  know  that  if  ds  =  0  in  the  K  reference  frame,  then  ds '  =0. 
Therefore,  in  general  ds  and  ds '  must  be  proportional,  i.e.: 

ds  =  a  ds '  (5) 
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where  a  can  depend  on  V,  i.e.,  a  =  a(V).  But  since  the  role  of  K 
and  K*  can  be  reversed,  i.e.,  K  is  moving  the  K'  is  not,  then 

ds ’  =  a  ds  (6) 

and  therefore  a  =  1.  Landau  and  Lifshitz  [Ref.  2]  give  a  more 
general  discussion  along  these  lines  including  among  other  things 
the  possibility  that  a  *  -1.  In  any  case  we  must  have  ds  =  ds ' 
and  therefore  s  =  s'.  Therefore,  the  interval  between  two  ents 

(s  or  ds)  is  the  same  in  all  inertial  reference  frames,  i.  it 
is  invariant  under  transformation  from  one  inertial  syster  o 
another.  In  particular,  this  invariance  is  the  mathematic 
expression  of  the  constancy  of  the  velocity  of  light.  Finally  we 
have : 

ds2  =  c2dt2  -  dx2  -  dy2  -  dz2 

=  c2 (dt ' ) 2  -  (dx  ' ) 2  -  (dz')2 

-  (ds')2  .  (7) 

To  obtain  the  mathematical  form  of  the  GLT  we  first  consider 
a  clock  fixed  in  the  K'  system  (dx '  *  dy ’  =  dz '  =0)  and 
therefore  moving  with  respect  to  the  K  system  along  the  x  axis 
with  velocity  V.  Now  this  clock,  viewed  in  the  K  system  would  be 
at  position  dx,  dy,  dz,  after  time  dt.  The  time  dt  would  be 
measured  with  an  identical  clock.  We  have  then: 


This  represents  the  "twin  paradox,"  i.e.,  moving  clocks  "appear 
to  run  slower  than  stationary  clocks. 


We  now  seek  the  transformation  of  an  event  recorded  in  the  K 
system  to  the  same  event  as  it  would  be  recorded  in  the  K' 
system.  We  restrict  the  discussion  to  motion  of  the  K’  system 
along  the  x  axis  of  the  K  system  as  one  can  always  rotate  the  x, 
y,  z  coordinates  to  a  system  with  the  motion  of  the  K'  system 
along  the  x  axis  of  the  K  system.  In  this  case,  the  y  and  z 
coordinates  are  unaffected  and  we  have  y  =  y’  and  z  =  z'.  Now  we 
can  seek  a  linear  transformation: 


"tT  fa  1 1 1 '  +  a12x' 
_xj  _a2  ict  '  +  a22x  - 


(ID 


Clearly,  we  must  have: 
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where  I  is  the  unit  matrix.  From  the  moving  clock  experiment 
above  we  have  immediately: 


a  1 1 


1 _ 

V1  -  (V/c)2 


(13) 


Now  consider  the  motion  of  the  K'  origin  as  seen  in  the  K 
system.  We  have  x'  =  y'  =  z'  =0.  and 


t  =  a  1 1 1 '  =  yt  ' 


x  =  a2 i ct  ' 


an  - 


1  x  =  V 
c  t  'c 


=  y6 


(14) 


Since  the  roles  of  K  a  K '  can  be  reversed,  we  must  have: 


I«ij(V)]  1  =  CaiJ (-V) ]  , 


(15) 


y  a12(-e) 

y  a l2 ( B  ) 

-Byc  a22(-6) 

Bye  a22 ( 3 ) 

(16) 
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which  immediately  gives: 


«12(*) 


=  JLl 


c 


a22  =  y 


Therefore  the  Lorentz  transformation  can  be  written: 


or 


Y  0Y/C 

t  ' 

C0Y  Y 

-  « 

X  ' 

(17) 


(18) 


V/e  can  write  these  expressions  in  a  more  symmetrical  form  by 
introducing  a  notation  that  considerably  simplifies  the  writing 
of  equations  and  makes  quite  clear  the  nature  of  the 
transformation.  We  deal  with  the  three  space  variables  as  a 
(contravariant )  vector: 


(19) 


(Since  we  are  dealing  with  flat  space-time,  we  need  not  be 
concerned  with  the  distinction  between  covariant  and 


contravariant  vectors  and  tensors.)  The  time  coordinate  can  be 
given  in  the  same  units  as: 


x 


0 


ct 


(20) 


Therefore  the  Lorentz  transformation  to  the  K'  system  moving  with 
velocity  V  along  the  x  axis  is: 


.0  , 


.1  , 


A  - 


,  3  , 


Y 

-By 

o 

o 


-By  o  0 

Y  0  0 

0  10 

0  0  1 


(21) 


In  the  following  we  will  label  the  spatial  (three-vector)  as  x, 
and  a  Lorentz  four-vector  as  'x"l  Therefore  the  Lorentz 
transormation  along  the  x  axis  is: 


x’  =  Li  (6)  'x 


(22) 


where 


Ll(fS) 


Y  -By  0  0 

-By  y  oo 

0  0  10 

0  0  0  1 


(23) 
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r  t. 


In  general  the  velocity  of  K'  will  be  in  an  arbitrary 
direction.  Therefore,  we  must  make  a  transformation  by  rotation 
of  the  three  spatial  axes  such  that  the  motion  of  K'  is  along  the 
positive  x  axis.  After  the  Lorentz  transformation,  we  can  then 

A 

reverse  the  rotation  to  obtain  the  general  transformation.  Let  U 
be  the  unit  three-vector  in  the  direction  of  motion  seen  in  the  K 
system.  We  can  chose  arbitrarily  the  other  two  unit  three-vector 
bases  for  the  transformed  system.  For  our  case  with  satellite 
trajectories  we  take  the  natural  base  vectors  as  the  across  track 

A  A 

(W)  and  approximate  radial  (R)  directions.  If  the  satellite 
position  and  velocity  three-vectors  are  r  and  r,  then: 


r«i 


u 


u2 

.uj 


l  rz — r 

r/y  r  •  r 


W  = 


"I 

W2 

LW3 


=  r  x 


-  J  -  -  -  - 

r/y  (r  x  r)  •  (r  x  r) 


R  = 


R1 

R2 

-R3J 


=  U  x  W 


(24) 


Then  the  desired  transformation  of  a  four-vector  is: 


where 


A(r,  r) 


10  0  0 
0  Ux  U2  U3 
0  Wj  W 2  W3 
0  Rj  R2  R3. 


and  the  inverse  is : 

10  0  0 

.  1  0  U,  »!  Rt 

A~(r ,  r)  = 

0  U2  W2  R2 
0  U3  W3  R3 


(2 


(2 


We  are  now  in  a  position  to  write  down  the  general  Lorentz 
transformation  to  a  moving  frame,  but  a  frame  with  its  space 
basis  three-vectors  parallel  to  the  original  frame.  This  is 
called  boost  without  rotation  by  Moller  [Ref.  4],  and  is 
analagous  to  Fermi-Walker  transport  in  the  General  Theory 
[Ref.  5],  This  transformation  is  simply: 


where  02  =  r  •  r/c2. 

With  a  little  algebra  one  can  write: 


-BrUi 


-eru2 


-0yU; 


S  =  S(U,  6)  = 


-BYUj  l+CY-lJUi2  (Y-llUiUz  (Y-lJUiUa 

-6yU2  (Y-l)UiU2  1+(y-1)U22  (y-1)U2U3 
-ByU3  (y-1)U!U3  (Y-1)D2U3  1+(Y-1)U32 


Note  that  this  transformation  is  in  fact  independent  of  the  basis 

A  A 

vectors  W  and  R,  as  it  should  be. 

We  now  define  the  frequency  four-vector  that  describes 
propagation  of  electromagnetic  plane  waves.  A  plane  wave 

A 

propagation  in  the  direction  of  the  unit  three-vector  i  is  a 
function  of  : 


r  •  l 


The  ti  are  the  direction  of  cosines  of  the  normal  to  the  plane 
wave  in  the  direction  of  propagation.  The  vector  potential  of 
such  a  wave  can  be  written  as  the  real  part  of  a  complex 
expression : 


A  =  Re{A0  exp(-ioi(t  - 


=  Re{A0  exp(-ij~(x°  -  x1*!  -  x2i2  -  x3*3))} 


which  can  be  written  in  terms  of  a  frequency  four-vector: 


as : 

A  =  Re  A0  exp  (i  4») 


where  the  phase  is: 


— 1 

l-» 

o 

o 

o 

_ 1 

~x0' 

0-100 

x1 

0  0-10 

x2 

o 

o 

1 

H* 

l _ 

_x3 

The  matrix 


[g  ] 

ij 


10  0  0 
0-100 
0  0-10 
0  0  0  -1 


(32) 


(33) 


(34) 


(35) 


is  in  fact  the  flat  space-time  metric,  and  the  occurrence  of 
gij  in  (32)  has  a  fundamental  significance  that  is  beyond  the 


14 


scope  of  this  discussion.  Now  in  a  moving  reference  frame,  the 
phase  would  be  viewed  as: 


♦  "  ♦*  “I’*1  *2*  I«ijl  * 

-lilt1  t2  t3]  [g±J  ]  T  (36) 

which,  using  (28),  can  be  written  as: 

[1  i1  t2  t3lg1J[S-1l  x7  (37) 

Now  using  the  relation: 

[gijltS]-1  =  [S][glJ]  (38) 

we  have 

♦  -  *  |  [1  l1  l2  t3][S]gij  P 

-  |'[1  t1'  l2'  t3']fgiJJ  F 

-  [V° '  V1 '  v2 '  v3 •  1  [g±J  ]  x7  (39) 


This  shows  that  the  frequency  four-vector  (32)  is  indeed  a 
four-vector  and  transforms  like  any  four-vector. 


The  frequency  four-vector  contains  all  the  information  about 
the  frequency  and  direction  of  propagation  of  a  plane  wave.  The 


Lorentz  transformation  describes  how  this  frequency  four-vector 
is  viewed  in  various  reference  frames.  Before  proceeding  with 
the  radar  doppler  we  can  now  easily  demonstrate  some  classical 
results,  which  can  provide  some  insight  into  the  deeper  meaning 
of  the  Lorentz  transformation. 

If  we  consider  the  geometrically  simple  situation,  where  K' 
system  is  moving  along  the  x  axis  of  the  K  system.  If  the 
direction  of  transmission  in  the  K  system  is  6,  then  the 
four-vector  in  the  K  system  is: 


r— i 
V  = 


v 


1 

cos  e 
sin  6 
0 


(40) 


This  four-vector  transformed  to  the  K'  system  is: 


- 

1  -  B  cos  9 

1 

-0  +  cos  0 

=  v ' 

cos  0  ' 

sin  0/y 

sin  0' 

o 

0 

The  apparent  frequency  is : 


v'  *  y v ( 1  -  e  cos  0)  (42) 


1C 


with  the  wave  moving  in  the  direction  9'  given  by: 


tan  9 ' 


v1  sin  9* 
v*  cos  9' 


sin  9 _ 

y(-B  +  cos  9) 


(43) 


Of  course  the  moving  frame  sees  the  wave  coming  from  the 
direction  0a'  =  n  -  0'.  The  moving  observer  would  deduce  that 
the  wave  was  actually  coming  from  the  direction  0a  =  it  -  9 
using : 


tan  9„  ' 
a 


sin  0_ 

_ a _ 

Y(6  +  cos  9  ) 

a 


(44) 


These  are  the  classical  expressions  for  the  doppler  shift  and 
aberration  as  seen  in  a  moving  frame. 

Now  we  sit  in  the  moving  (K')  system,  and  reflect  the 
incoming  wave  ( v ' )  back  along  the  direction  of  propagation 
(rvJ.).  This  could  be  done  with  a  cube  corner.  Mathematically 
this  is  expressed  as: 


If  we  then  view  this  reflected  wave  in  (i.e.,  transform  to)  the 
stationary  (K)  system  we  obtain: 
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I 


1-26  cos  9  +  62 

_  „  -cos  9(1  +  B2)  +  26 

V  =  LjC-B)^*  =  vy2 
r  r  -sin  0/y2 

0 


(46) 


Now  the  frequency  is  perceived  to  be  shifted  by: 


v  -  v 

- -  =  y2  ( -23  cos  0  +  232  ) 


and  propagated  in  the  direction  0r: 


(47) 


tan  0 

r 


_ sin  0 _ 

Y 2 ( cos  0(1  +  B2 )  -  23) 


appearing  to  come  from  the  direction: 


sin  0 

tan  0  =  - 

r  y2(cos  0a(l  +  32 )  +  23) 


(48) 


(49) 


Note  that  this  is  not  the  original  direction  of  outward 
propagation  (+180  deg)  in  the  K  system,  though  the  wave  was 
reflected  along  the  transmitted  direction  in  the  K'  system.  This 
effect  is  known  as  velocity  aberration.  In  this  case  the 
refected  wave  will  not  return  to  the  transmitter.  For  a 
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satellite  at  the  height  of  Lageos  (7503901,  8820)  the  return 
would  be  120  m  away  from  the  transmitter.  Therefore  the  cube 
corners  mounted  on  Lageos  and  other  satellites  for  laser  ranging 
must  be  imperfect,  spreading  the  return  pulse  such  that  some 
energy  will  return  to  the  transmitter  site. 

From  the  point  of  view  of  the  moving  frame  (K')  one  would 
have  to  reflect  the  plane  wave  along  the  direction: 


-  8  cos  0 ) 


-(cos  0  +  8)/(  1+8  cos  9) 
-sin  0/( y  [  1  +  8  cos  0]) 


to  return  to  the  transmitter.  In  this  case  the  apparent 
frequency  would  be: 


V  =  vy2(l  -  28  cos  0  +  282  cos2  0) 
r 


One  could  also  return  the  plane  wave  along  the  direction: 


v'  *  yv(l  -  8  cos  9) 


-(8  +  cos  0(1  -  26z))/(l  -  8  cos  6) 
-sin  9/y 


F 


which  would  exactly  cancel  the  effects  of  Special  Relativity  in 
observed  frequency.  One  notes  that  all  these  effects  are  of 
second  order  in  the  small  quantity  B.  This  is  why  the  effects  of 
TSR  are  experimentally  difficult  to  observe,  and  that  some  care 
must  be  used  in  the  details  of  computation. 

These  examples  assume  the  moving  target  controls  the 
direction  of  reflection;  more  precisely  the  direction  of 
retransmission.  In  the  case  of  radar  reflection,  the  target  is 
expected  to  return  some  radiation  to  all  directions  in  the  half 
space  normal  to  the  incoming  wave  direction  in  the  stationary 
system.  One  can  express  the  returned  four-vector  in  terms  of  the 
four-vector  in  the  stationary  system.  That  is  how  (50)  was 
obtained . 

Within  this  framework  we  can  now  formulate  the  radar  Doppler 
equation.  Recall  that  we  ultimately  need  to  write  expressions 
for  the  transmitted  and  observed  frequency  in  terms  of  the 
inertial  geocentric  position  and  velocity.  Let  us  introduce  the 
following  notation.  The  time  of  transmission,  the  time  of 
reflection  at  the  target,  and  the  time  of  reception  are 
respectively:  tl,  t2,  t3  as  defined  by  a  clock  in  our  inertial 

geocentric  system.  Prom  (10)  we  see  that  the  (TSR)  correction  of 
the  time  measured  at  the  transmitter  and  receiver  to  inertial 
geocentric  time  is  negligible.  The  position  and  velocity  of  the 
transmitter  at  tl  are  denoted  are  Ri,  ,  the  position  and 
velocity  of  the  target  at  t2  are  r2 ,  r2,  and  the  position  and 
velocity  of  the  receiver  are  R3,  R3 ;  all  expressed  in  the 
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inertial  geocentric  reference  frame.  The  unit  three-vector  of 
each  velocity  is  denoted  U.  the  magnitude  of  each  velocity 

vector  is  8  in  units  of  the  speed  of  light  (c).  Therefore  we 

•  •  • 

_  A  _  A  _  A 

have  Ri  =  c8iUlf  r2  =  082^2  and  R3  =  C83U3.  We  have  the  vector 
from  transmitter  to  target  Pi2  =  r2  -  Ri  and  from  target  to 
receiver  p2 3  =  R3  -  r2 .  The  ranges  are  (pi2)2  =  Pi2  •  Pi2  and 
(P23)2  =  P23  *  P 2 3  *  Now  we  can  define  the  four-vector  wave 
propagating  between  transmitter  and  target  in  the  inertial 
reference  frame  as : 

where 

A  _ 

*12  ~  P12/P12 

Now  this  wave  as  seen  by  the  (moving)  transmitter  is: 


=  S(Uj  ,  6 1)1^ 

and  the  target  sees  this  wave  as: 


(54) 


i- 

L _ £ 


=  S(U2,  62)^ 
Using  (54)  we  have: 
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(55) 


•ii.-...,  ■-  aatdfciftE 


(56) 


v"1  =  S(U2,  e2)S(ult  -0i)^ 

The  same  reasoning  relates  the  received  frequency  four-vector  as: 

vj  -  s(u2,  e2)S(u3.  -e3)vj  (57) 

Since  the  transmitted  frequency  seen  hy  the  target  (v^°")  and 
the  reflected  frequency  at  the  target  (vr°")  are  the  same,  we 
can  equate  the  first  element  of  these  four-vectors  and  obtain  the 
desired  relation  between  the  transmitted  (v't)  and  received 
(v'r)  frequency  in  their  respective  frames. 

For  calculation  we  need  not  compute  the  whole  Lorentz 
transformation  matrix.  We  can  write  out  the  elements  of 


S(U,  0)S(U',  6')  =  [S.j] 

as : 


(58) 


oo 


Soi 


io 


ii 


A  A 

YY  '  +  0  0 1 YY ' ( b  *  U  '  ) 

-Y[0'Y'Uj  +  01L  +  a'0Ui(U  •  U')] 

-Y,[0YUi  +  0'U^  +  a0  '  IL  (U  *  U '  )  ] 

1  +  «U2  +  a’(U^2)  +  UiU^(0Y0,Y'  +  oa ’ (U  •  U’)] 

00'YY'U.U'  +  aU.U.  +  a'U'.U'  +  aa '  U !  U,  ( U  •  U') 

J  ^  J  ^  J  ^  J  i 

BB'YY’U  O'.  +  aU.U.  +  a  '  U  ’  U  '.  +  aa’U'U.(U  •  U’) 

^  J  ^  J  ^  J  ^  J 


(59) 

i  =  1,2,3 
i  =  1,2,3 
i  =  1,2,3 

i  >  j  >/  1 
j  >  i  >,  1 
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where  Ui  are  the  three  components  of  V  and  U'i  are  the  three 
A  2 

components  of  U',  and  a  =  (y  -  1)  *  3  /2  is  a  small  quantity  of 
second  order.  Note  that  we  have  used  positive  0'  in  (58). 
Therefore  in  (56)  6'  =  -3i,  and  in  (57)  6'  =  -63. 

.  A  A 

Let  us  call  o  i j  =  Sij(U2,  32 >  ul »  -&i)  and 
S3jj  =  Sjj(U2,  32  ,  U3  ,  -33).  Then  we  have: 


v ' 
r 


[S3.  ] 
1  01  1 


[S1 
1  00 


3 

I 

i=l 


01 


*12  1 


(60) 


or : 


00 


*  £  sii 


v  '  =  v 


i=l 


t  3 

s3  +  y  s3. 

00  i  =  l  01 


(61) 


The  frequency  shift  seen  is  therefore: 


v'  -  v' 
r  t 


00 


s3  +  y  s1.  ti2  -  y  s3  1 2 3 

00  i=l  oi  12  i=l  oi  23 


00 


3 

+  l 

i=l 


,3 

°Oi 


*  23 


(62) 


This  expression  is  the  desired  result,  without  approximation,  and 
no  assumption  that  the  transmitting  and  receiving  site  are  the 
same.  For  purposes  of  computation  (62)  should  be  used.  To  see 


the  meaning  of  the  terms,  we  note  the  following.  First,  B'  is 
the  velocity  of  the  transmitter  (receiver  or  station)  and  is  in 
general  smaller  than  0  for  ground  based  observations.  If  we 


consider  o'  ■  8,2/2  negligible  and  if  we  therefore  take  o'  *  1, 
then : 

soi  =  ^ 0  *  ui  +  1  =  1’2>3  <63> 

Recall  that  S'  =  8i  or  -B3  and  we  see  that  the  Lorentz 
transformation  depends  on  the  difference  in  velocity  of  the 
observer  and  target,  which  is  consistent  with  the  first  tenent  of 
TSR.  Further,  if  we  assume  the  station  motion  is  rectilinear 
(i.e.,  unaccelerated)  then  Bx  =  83,  Uj  =  U3  and  S^q  =  S300 
and  S'0i  =  S30i.  Therefore  the  observed  frequency  shift 
becomes : 


3 

1 

i-1 


oi 


<*I2  ~ 


*2  3  ) 


OO 


3 

I 

i=l 


S'  . 
Oi 


„i 
l2  3 


(64) 


This  approximation  ignores  the  Lorentz  boost  due  to  the  station 
velocity,  but  preserves  the  different  pathes  of  transmission.  If 
finally  one  takes  the  mean  station-target  vector 


*22  =  (^2  -  ^2  ) / P 2 2 


(65) 
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k. 


j  -a  11 '  * 


with 


*12  ”  *22  “  -*23 


(66) 


and 


Soo  =  Y’  Soi  =  "Y(r  ~  R)/c 


(67) 


we  have : 


v  '  -  v ' 
r  t 


-2(r  -  R)  •  l22 


c  ( 1  + 


(r  -  K) 


l22 


yc 


(68) 


or 


c 

2 


v  '  -  v. 


=  “(r  -  R)  *  t22  +  t(r  -  H)  •  *22]2/c 


(60) 


the  classical  approximation. 

If  one  wants  an  approximation  with  the  accuracy  of  (64), 
but  in  terms  of  the  mean  geometric  station-target  vector  one  has 


v  '  -  v  1 
c  r  t 


■  -(r  -  R)  •  *22  +  t(r  -  H)  •  * 22 1 2 /c 


-(r  -  R)  •  AR/ p 


(70) 


where  AR  =  R2  -  Ri  =  R3  -  R2. 

Since  the  second  and  third  terms  of  (70)  are  terms  of  second 
order  they  could  be  calculated  with  a  predicted  trajectory.  We 


can  then  define  the  mean  geometric  range  rate  in  terms  of  the 
observable  as: 


*22 ]2 /c 


+  (r  -  H)  •  £K/p 


or 


+  (r  -  1?)  •  TR/p 


(71) 


(72) 


HI.  DISCUSSION 

In  all  cases  we  have  expressed  the  frequency  shift  as  viewed 
in  a  moving  frame  (the  observable)  in  terms  of  positions  in  our 
interial  geocentric  reference  frame  (the  computable).  Therefore, 
the  vector  quantities  in  (62),  (64),  and  (69)  are  considered 
geometric  quantities  and  are  the  three-space  quantities  used  with 
Newtonian  Mechanics.  The  correction  terms  for  the  Newtonian 
three-space  are  small.  The  elaboration  of  a  trajectory  is 
considerably  simplified  if  done  in  three-space,  allowances  being 
made  with  small  corrections.  Therefore,  the  meaning  of  (62), 
(64),  and  (69)  is  as  follows.  A  trajectory  or  orbit 
determination  is  made  in  convenient  Newtonian  three-space.  Radar 
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observations  are  used  to  determine  the  orbit  state  vector 


iteratively  by,  say,  the  method  of  Least  Squares.  In  this 

calculation,  the  observed-computed  difference  is  needed.  In  this 

case  (62),  (64)  or  (69)  is  used  to  obtain  the  computed  frequency 

shift.  Then  approximating  the  difference  as  an  error  in  the 

•  • 

A 

geometric  range  rate  (r  -  H)  •  122  an  observation  equation  is 
formed  to  improve  the  satellite  state  vector. 

The  discussion  thus  far  has  employed  the  TSR.  We  know  that 
the  space-time  is  in  fact  not  flat  but  curved  due  to  the  presence 
of  mass  or  energy.  This  is  known  as  the  General  Theory  of 
Relativity  (GTR).  For  our  purpose  we  can  approximate  the  local 
space-time  curvature  with  the  Schwarzschild  Metric  [Ref.  5], 

There  will  be  effects  due  to  the  Sun  and  the  Earth.  We  know  that 
based  on  the  GTR,  there  are  a  number  of  effects  that  must  be 
considered.  First,  the  dynamics  of  a  particle  must  be  modified. 
This  is  most  easily  accomplished  with  a  small  "correction"  term 
to  the  Newtonian  model.  For  close  earth  satellites,  this  is 
trivially  accomplished  [Ref.  6].  Second,  the  velocity  of  light 
is  slower  in  the  presence  of  a  gravity  field,  than  it  is  in  flat 
space  time.  Using  a  light  travel  time  of  electromagnetic 
propagation  to  measure  distance  therefore  requires  a  correction. 
In  this  case  the  effect  of  the  earth's  mass  is  very  small,  but 
the  sun's  mass  is  not.  The  correction  can  be  estimated  as 
follows.  The  space-time  relation  between  the  proper  time  of 
radial  transit  of  a  photon  (dt)  and  the  metric  distance  (r  -  r0) 
Ln  the  Schwarzschild  metric  can  be  shown  to  be 
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(73) 


r 


cdt  =  r  -  r0  +  2— £  ln(r/r0) 
c 

where  ro  is  the  radius  of  the  transmitter,  r  is  the  radius  of  the 
target,  and  GM  is  the  gravitational  constant  multiplied  by  the 
mass.  GM/c2  has  the  units  of  length  and  is  sometimes  referred  to 
as  the  Schwarzschild  radius.  The  first  two  terms  of  the  right 
hand  side  of  (73)  are  the  flat  space-time  (Newtonian)  distance, 
and  the  last  term  can  be  viewed  as  a  General  Relativity  term. 
Examples  of  the  correction  are  given  in  Table  I. 

TABLE  I 

SPACE  TIME  CURVATURE  EFFECT  ON  RANGE  MEASUREMENT 


Mass 

Range  Measurement  To 

Correction 

Earth 

Lageos 

0.6 

cm 

Earth 

Synchronous  Satellite 

1.6 

cm 

Earth 

Moon 

3.6 

cm 

Earth 

Sun 

8.9 

cm 

Sun 

Lageos 

12.5 

cm 

Sun 

Synchronous  Satellite 

62.8 

cm 

Sun 

Moon 

741.0 

cm 

Sun 

Sun  (lf^!)) 

16.9 

km 

Now  the  conventional  value  of  the  speed  of  light  in  vacuum  (c  = 
2.99792458  *  0.000000012  x  10 10  cm/sec)  is  based  on  laboratory 
measurements  at  the  earth's  surface.  Therefore  measurements  of 
distance  near  the  earth  implicitly  use  the  appropriate  value  of 
c.  in  other  words,  a  measurement  of  dt  at  the  earth's  surface 
correspond  to  r  -  ro  if  the  absolute  speed  of  light  is  reduced  by 

i£  a  2-SM_ 
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Therefore,  the  speed  of  light  at  infinity  is  greater  by  ths 


value.  For  M©  and  a®,  6c/c  =  2  x  10“8.  We  can  summarize  the 
situation  as  follows.  The  choice  of  a  particular  value  of  c 
depends  on  the  reference  frame  chosen.  Discussion  of  data  in  the 

near  earth  environment,  i.e.,  where  one  can  ignore  the  change  in  i 

i 

c  due  to  the  change  in  distance  from  the  sun,  the  constant  ' 

i 

i 

obtained  in  terrestial  laboratory  measurements  is  correct.  In  j 

this  case,  only  the  variation  in  the  space  time  metric 
represented  by  the  earth's  mass  should  be  used.  Since  ranges  to 

I 

the  moon  need  be  corrected  by  less  than  5  cm,  this  can  be  ignored  ! 

in  most  applications.  However,  when  considering  general  solar  I 

I 

system  measurements,  then  a  frame  with  c  defined  at  infinity  must 
be  used.  The  metric  distances  in  the  earth's  environment  will  of 
course  then  have  to  be  expressed  in  this  new  system.  Therefore, 
for  earth  orbiting  bodies,  within  the  moons  orbit,  the 
corrections  given  in  Table  I  are  not  applicable.  Finally,  the 
frequency  of  our  standard  clocks  will  change  in  a  potential 
field.  This  is  potentially  the  most  significant  effect  with 
Doppler  measurements.  The  appropriate  description  includes  a 
transformation  from  station  to  satellite.  This  transformation  is 
the  change  in  frequency  along  the  geodesic  path  from  station  to 
satellite. 

It  is  found  that  the  frequency  would  change  by  a  factor 
1  +  A$,  where  A$  is  the  difference  in  potential  between  the 
transmitter  and  target.  The  reverse  change  takes  place  on  the 
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return  path.  If  the  transmitter  and  receiver  are  at  the  same 
potential  (e.g.,  at  the  same  geopotential  height)  then  the  change 
up  and  down  would  be  exactly  the  same.  However,  if  the 
measurement  is  of  an  oscillator  onboard  the  target,  then  such  a 
term  should  be  included  in  (62)  [Ref.  8].  The  consequences  of 
this  effect  for  multistation  observatories  must  be  examined 
further. 
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